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ABSTRACT 

Context. HH 223 is a knotty, wiggling nebular emission of ~ 30" length found in the L723 star-forming region. It lies projected onto 
the largest blueshifted lobe of the cuadrupolar CO outflow powered by a low-mass YSO system embedded in the core of L723. 
Aims. We analysed the physical conditions and kinematics along HH 223 with the aim of disentangling whether the emission arises 
from shock-excited, supersonic gas characteristic of a stellar jet, or is only tracing the wall cavity excavated by the CO outflow. 
Methods. We performed long-slit optical spectroscopy along HH 223, crossing all the bright knots (A to E) and part of the low- 
brightness emission nebula (F filament). One spectrum of each knot, suitable to characterize the nature of its emission, was obtained. 
The physical conditions and the radial velocity of the HH 223 emission along the slits were also sampled at smaller scale (01' 6) than 
the knot sizes. 

Results. The spectra of all the HH 223 knots appear as those of the intermediate/high excitation Herbig-Haro objects. The emission 
is supersonic, with blueshifted peak velocities ranging from -60 to -130 km s"'. Reliable variations in the kinematics and physical 
conditions at smaller scale that the knot sizes are also found. 

Conclusions. The properties of the HH 223 emission derived from the spectroscopy confirm the HH nature of the object, the super- 
sonic optical outflow most probably also being powered by the YSOs embedded in the L723 core. 

Key words. ISM: jets and outflows — ISM: individual objects: L723, HH 223 — stars: formation 



1. Introduction 

Lynds 723 (L723) is an isolate dark cloud located at a distance 
of 300+150 pc (Goldsmith et al. 1984) that shows evidence of 
low-mass star formation, where one of the few known cases of 
_ a quadrup olar C O outflow (two separate pairs of red-blue lobes; 
k> , Lee et al. l2002l and references therein) has been report ed. Th e 
j_] ■ 3.6 cm radio continuum source VLA 2 (Anglada et al. Il996h . 
d ' towards the centre of the CO outlfow, seems to harbour the ex- 
citing outflow source. VLA 2 is embedded in high-density gas 
traced by NH3, showing evidence of gas heating and line broad- 
ening (Gkart et al. [1997i) . The dense envelope of the source has 
als o been observed at submil limete r wavelengths by Shirley et 
al. (I2OO2I) . and Estalella et al. (|2003|) . R ecent w orks by CaiTasco- 
Gonzalez et al. (11)08) and Girart et al. ( |2008|) . show that VLA 2 
is a multiple system. Carrasco-Gonzalez et al. ( 120081) detect at 
least four young stellar objects (YSOs) and propose that the mor- 
phology of the CO outflow is actually the result of the superpo- 
sition of three independent pairs of CO lobes. They also propose 
that one of the YSOs (VLA 2A) is exciting the largest pair of 
CO lobes and the system of emission-line nebulosities, reminis- 



cent o f Herbig-Haro (HH) objects, first reported by Vrba et al. 
(fT986l) . 



Send offprint requests to: R. Lopez 

* Based on observations made with the 2.6 m Nordic Optical 
Telescope operated at the Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias. 



In a previous work (Lopez et al. |2006!), we presented deep 
narrow-band Ha and [S 11] images of L723, which confirmed the 
emission-line nature of the nebulosities. The Ha image, with 
the CO outflow contours from Lee et al. (2002) superposed, is 
shown in panel (a) of Fig.[T] The higher angular resolution of the 
images allowed us to resolve the detailed structure of HH 223, 
the br ightest "linear emission feature" reported by Vrba et al. 
(m^. HH 223 appeared as a set of knots embedded in a fainter 
nebula, following a wiggling pattern, reminiscent of a stellar jet, 
that is probably part of a parsec-sc ale HH outflow as suggested 
by CO observations of Lee et al. (|2002|) . o ptical (Lopez et al. 
20061) and near-infrared (Palacios & Eiroa [19991 : Lopez et al. 
2009") imaging of the shocked emission. However, there were no 
kinematic data to establish the nature of the line emission, i. e. 
whether the emission is supersonic (as in typical stellar jets) or, 
in contrast, is nearly stationary as would be expected for a wall 
cavity. With the aim of exploring the kinematics and physical 
conditions in HH 223, we performed long-slit spectroscopy of 
the object covering all its bright knots (A to E) and a part of the 
low-brightness nebula (F filament). The results are presented in 
this paper. 
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Fig. 1. (a) Ha image of the L723 field. HH 223 and the other Ho- nebulosities given in Lopez et al. ( |2006|) have been labeled and 
enclosed within boxes. The triangles mark the positio ns of the H2 knots Kl and K2 of Palacios & Eiroa (Il999h . Contours of the 
high-velocity CO (7=1— >0) emission from Lee et al. |2002| have been superposed (the das hed co ntours correspond to redshifted 
emission and the solid contours to blueshifted emission). The VLA2 source (Anglada et al. Il996h is marked, (b) Close-up of the 
image showing HH 223. The slit p ositions for the spectra acquisition (SI, S2) are indicated. Knots are labeled according to the 
nomenclature of Lopez et al. (l2006h . 



2. Observations and Data Reduction 

Spectra of HH 223 were obtained with the 2.6-m Nordic 
Optical Telescope (NOT) of the Observatorio del Roque de los 
Muchachos (ORM, La Palma, Spain) on 13 July 2007. The 
ALFOSC (Andalucfa Faint Object Spectrograph and Camera), 
with an image scale of 0'.'188 pixel"', was used in its spec- 
troscopic mode. The 600R grism # 8, covering the wavelength 
range 5825-8350 A with a dispersion of 1 .3 A pix"' (spectral res- 
olution ^ 60 km s ') was used. In order to cover all the HH 223 
bright knots, the spectra were acquired through a long-slit of 6.'5 
length and 1" width, positioned at two position angles (see panel 
(b) of Fig.[B: at PA 73°, along knots A and B (SI); and at PA 
79°, along knots C to F (S2). Individual exposures of 1800 s 
each were made to get total exposure times of 5400 s and 7200 s 



for the SI and S2 slit positions, respectively. Special care was 
taken in achieving accuracy to point the slit on the knots. The 
setup only allows horizontal or vertical slit positions. Thus, a 
shorter exposure (200 s) of the field oriented at the desired PA 
was previously acquired through a narrow-band Ha filter. Then, 
an image of the slit on the detector was made and the proce- 
dure was iterated twice to achieve an accurate slit position be- 
fore acquiring the spectrum. Seeing was 0'.'6 for most of the run. 
The data were processed with the standard tasks for long-slit 
spectroscopy of the IRApQ package, which include bias sub- 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 2. Spectra of the HH 223 knots from east (A, top-left panel) to west (F filament, bottom-right panel) obtained by averaging the 
signal inside the slit crossing the knot. In all the spectra, the intensity has been normalized to the Ha peak intensity. The emission 
lines detected in each spectrum have been labeled in the panels. The slit position (and the size of the region averaged along the slit 
to get the spectrum) are indicated in each panel. 



straction, flat-fielding correction, wavelength calibration using 
spectra from a Ne lamp acquired at several elevations for each 
slit position, and sky substraction. The achieved accuracy for the 
wavelength calibration, checked from the sky lines, was better 
than ~ 0.2 A (~ 10 km s"'). The accuracy in the determination 
of the position of the line centroids is ~ 0.2 A (~ 10 km s"') for 
the strong emission lines observed. The final spectrum through 
each slit position was corrected for cosmic ray events by median 
filtering all the exposures, after wavelenght calibration, obtained 
at the given slit position. The spectra were not flux calibrated. 

3. Results and Discussion 

3.1. The spectra of the HH 223 knots: Physical conditions 
and kinematics 

We obtaine d the spectrum for each HH 223 knot defined in 
Lopez et al. ( 2006) by averaging the emission within the slit win- 
dow crossing the knot. The six spectra obtained in such a way 
are displayed in Fig.|2l As can be seen from the figure, the char- 
acteristic HH emission lines from [Oi] A 6300, 6364 A, [Nn] 
A 6548, 6583 A, Ha and [S ii] A 6716, 6731 A were detected in 
all the knot spectra. In addition, other emission lines redwards of 
7000 A were detected in several knots: the [Fe ii] 14F multiplet 



lines (A 7155, 7173, 7388, 7453 A) were detected in knot A; the 
brightest [Feii] A 7155 A component of 14F multiplet was also 
detected in knots B and C, and the [On] A 1321, 7331 A and 
[Can] /I 7291, 7324 A fines (the two doublets being blended) 
were detected in knots A to E. 

3.1 .1 . Density and excitation conditions of tine knots 

From the spectra obtained as mentioned above, we explored the 
electron density (We), from the [S ii] 6716/6731 line ratio, and the 
excitation conditions, from the [N ii]/Ha' and [S iil/Ha line ra- 
tios. The fluxes were uncorrected for reddening. However, these 
line ratios are insensitive to differential extinction because of the 
wavelength proximity of the lines. 

Figures [3] and |4] show the location of the HH 223 knots in 
current diagnostic diagrams from line ra tios of nebular emission 
lines (see e.g. Canto 1981; Raga et al. Il996h . As can be seen 
from these figures, the line ratios of all the HH 223 knots lie 
in the diagrams region where are located the HH objects having 
an intermediate- or high-excitation spectrum. The diagrams also 
show that there are appreciable differences in both, excitation 
and density, among the knots. Knot A has the highest excita- 
tion spectrum. Concerning the rest of the HH 223 knots, there 
is some trend of increasing excitation going towards the west. 
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Table 1. Vlsr line centroids' of the HH 223 knots 



Knot 


[Oi] 6300 




[Nil] 6583 


[Sn] 6731 


Average 


A 


-80 ± 20 


-108 + 15 


-83 ± 15 


-108 ± 15 


-95 ± 33 


B 


-115± 15 


-122 ± 15 


-131 ± 15 


-128 ± 15 


-124 + 30 


C 


-23 ± 15 


-76± 15 


-67 ± 15 


-57± 15 


-56 ± 30 


D 


-80± 15 


-117 ± 15 


-92 ± 15 


-106 ± 15 


-99 ± 30 


E 


-63 + 15 


-99 ± 15 


-79 ± 20 


-79 ± 25 


-80 ± 33 


F 


-24 ± 15 


-26 + 15 


-1 ± 15 


-17 + 25 


-17 + 35 


Notes: 



(1) Derived from the spectra integrated within the slit aperture sizes indicated in Fig.|2] 
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Fig. 3. Position of the HH 223 knots in a diagnostic diagram of 
Hq'/[S ii] vs. Hq'/[N ii] (line ratios, obtained from the knot spectra 
of Fig.|2]l. The regions of the Planetary Nebulae (PN), Supernova 
Remnants (SNR), H II regions (H II) and Herbig-Haro objects 
(HHs), adopted from the figures by Canto (1981.) . are outlined. 
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Fig. 4. Position of the HH 223 knots in the excitation vs. density 
diagnostic diagram. The long-dashed line gives the separation 
between hi gh/int ermediate and low excitation HH spectra (see 
Raga et al. ll996h . The regions of the Planetary Nebulae (PN), 
Supernova Remnants (SNR) and H II regions (H II), adopted 
from the figures of Meaburn & White (il98Z) . are outlined. 



Values for We of the knots were derived from the [S ii] 67 1 6/673 1 
ratio, using the TEMDEN task of the IRAF/STSDAS package, and 
assuming Tg = 10'* K. The results are shown in Fig.|4] The elec- 
tron density is significantly higher in the brighter, southern knots 
A and B than in the rest of the HH 223 knots, the difference be- 
tween the highest (in knot A) and the lowest (in knot C) «e values 
being about one order of magnitude. Knots D, E and F have sim- 
ilar He values. Note, however, that knots B and C have spectra of 
similar excitation, in spite of their derived being significantly 
difi'erent. 

3.1 .2. Radial Velocities of the knots 

The radial velocitie^ of the HH 223 knots were obtained from 
the line centroids of Gaussian fits to the emission lines of the 
spectra of Fig.|2l Thus, the derived velocities correspond to the 
knot emission integrated within the slit aperture and do not ac- 
count for differences in velocity at spatial scales smaller than the 
knot size. In Table [T| we hst the velocities obtained from several 
emission lines and the average velocity. The velocities derived 
for the knots appear blueshifted relative to the ambient gas, the 
gas intersected by the slit moving supersonically towards the ob- 



^ All the velocities in the paper are referred to the local standard of 
rest (LSR) frame. A Vlsr -+10 -9 km s"' for the parent cloud has been 
taken from Torrelles et al. 119861 



server. Thus, the HH 223 knot emission properties are those ex- 
pected from an optical outflow. Note that the velocity derived 
for HH 223-F is significant lower. However, the spectrum of 
HH 223-F mostly arises from gas of the low-brightness nebula 
in which the knots are engulfed (F filament) instead of from a 
truly, compact knot. 

More detailed maps are shown in Figs. |5] and |6] For knot 
A, the position-wavelength (PV) maps of Fig.|5]show a velocity 
gradient within this knot, with higher blueshifted velocities to- 
wards the east, the trend being found in all the lines mapped. In 
fact, the PV map of the Ho' line suggests that there are two veloc- 
ity peaks, at Vlsr- -145 and -85 km s"', with a position offset 
^ 0'.'5. Interestingly, the naiTow-band Ha image also shows two 
intensity maxima within knot A (labeled Al and A2 in Fig. [T] 
upper panel), the positions of these intensity maxima being in 
good concordance with the positions of the two velocity compo- 
nents of the PV maps. This suggests that the emission labeled as 
knot A, extending over ~ 3'.'5 along the slit position SI, would 
be in fact a more complex structure composed of two emitting 
clumps, partially resolved by our data, each clump moving at a 
different blueshifted velocity. This is better seen in Fig.|7] where 
we show the PV maps for different lines (lower panels) and the 
intensity contours obtained from the Ha image, oriented accord- 
ing to the slit position S 1 . Note in addition that the peak velocity 
position in the Ha line (i. e. around the A2 substructure) is dis- 
placed from the peak velocity in the [Sii] A 6731 A line (i.e. 
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Fig. 5. Wavelength-position diagram of the HH 223 emission 
along slit position SI (PA 73°) for the lines indicated in each 
panel. The rest wavelength that coiTesponds to Vlsr^ km s ' 
is marked with a dashed line. On the [O ii] + [Ca ii] panel, the rest 
wavelengths marked correspond to the [Can] A 7291 A (dashed) 
and to the [On] A 7321 A (dotted-dashed) component of their 
respective doublets. The reference offset position corresponds to 
the [S ii] knot A peak intensity. The knots identification is labeled 
in the [O i] panel. 
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Fig. 6. Same as Fig.|5l but for the slit position S2 (PA 79°). The 
reference offset position corresponds to the [S ii] knot C peak 
intensity. The knots identification is labeled in the Hot panel. 



around the Al substructure). Spatial displacements between the 
Ha and [S ii] peak positio ns are also found in other stellar jets 
(e.g HH 110, Riera et al. 120031: HH 119, Galfalk & Olofsson 
12007 *). in which the supersonic jet outflow propagates in an in- 
homogeneous, quite dense ambient environment. 

3.2. Small-scale properties of the HH 223 emission 

Our data, having better spatial resolution than the knot sizes 
(ranging from ~ 2" to 3'.'5), suggest that there are changes in 
both kinematics and physical conditions at a smaller spatial 
scale. With the aim of looking for such variations, we extracted 
a set of spectra by co-adding the signal of each three adjacent 
pixels (i. e. by doing a spatial binning of ~ 0'.'6, of the order of 
the seeing) to cover without gaps all the HH 223 emission inter- 
sected by the slits. The spectra, still keeping a very good signal- 
to-noise ratio (several of them are displayed in Fig. [8]i are then 
useful for sampling the gas conditions with a step of ~ 0'.'6. 

The spectra displayed in the two upper rows of Fig. |8] cor- 
respond to the emission at adjacent positions through knot A 
(from east to west). A quick look to the panels shows apreciable 
changes in the relative intensity of the emission lines, suggestive 
of changes in excitation and electron density even by only mov- 
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Fig. 7. Upper panel: Contour plot of the Ha emission in knots 
A and B obtained from a narrow-band image of the L 723 field. 
Lower panel: PV maps for different emission lines, labeled in 
each panel. Contours are in counts (image and spectra are not 
flux calibrated) and the contour spacing is 10 per cent of the 
knot A peak intensity. All the panels are oriented at a PA=73° 
(corresponding to the S 1 slit position, crossing the HH 223 knots 
A and B in the east-west direction). The reference offset position 
is as in Fig. |5] The dashed lines mark the peak positions in the 
PV maps of the two substructures, Al and A2, of knot A. 



ing ~ 0'.'6 along the slit. Fig. [8] also displays some of the spectra 
obtained at several selected positions along slits S 1 (within knot 
B) and S2, for positions belonging to knots C, D and the F fila- 
ment. It should also be noted that the [S ii] /1 67 16 A line appears 
brighter than the A 6731 A line in the spectra extracted along 
S2, while the opposite occurs along S 1 . Given that the [S ii] A 
6716/6731 line ratio is a tracer of the density of the ionized gas, 
this finding is indicative of a general trend, both from the knots 
and from the interknot low-brightness regions, of the gas being 
less dense towards the northwest, i. e. towards the optical ob- 
scured region. 

Let us now examine the small-scale changes in the physical 
conditions using the line ratios derived from these spectra. 

3.2. L Small-scale excitation and density along HH 223 

Figures |9] and [TO] (lower panels) display the spatial behaviour 
along the slit for several line ratios, tracing the excitation and 



ionizaton conditions, and for the electron density. To help with 
the knot identification, the Ha image is displayed in the upper 
panel of these figures. 

As can be seen from Fig. |9] all the line ratios differ from 
knot to knot, but also vary at a smaller spatial scale, within the 
knot. The gas excitation is traced by the [N ii]/Ha and [S ii]/Ha 
line ratios. All the values derived along knots A and B corre - 
spond to the intermediate/high excitation HHs (Raga et al. 19961). 
However, the excitation is clearly different in both knots: the 
[S n]/Ha line ratios along knot A are < 0.6, while they are > 
0.9 along knot B. Note in addition that the lowest [S ii]/Hq' line 
ratio value (thus the highest excitation) is found inside knot A, 
at positions offset > +1" from the two knot A peaks. In contrast, 
the [S iil/Ha line ratio values of knot B are the highest derived 
for all the sampled emission, thus indicating that the lowest ex- 
citation in the emission is found around knot B. Th ese res ults 
are fully consistent with those derived by Lopez et al. ( l2006l) us- 
ing the [S ii]/Ha ratio map obtained from the narrow-band CCD 
images. 

The behaviour of the electron density through knot A and B, 
derived from the [S ii] A 6716/6731 line ratio (see Section 3.1.1), 
is also shown in Fig.|9l An electron density «e - 200-300 cm"^ 
is found at positions away from the knot peaks, coinciding with 
the regions where the slit mainly intersects gas coming from the 
low-brightness emission. In contrast, values in the range 1200 < 
We < 3000 cm"^^ are found through knot B and 1500 < < 2500 
cm"^ through knot A, except for the position around the A2 sub- 
structure, where the derived [Sii] 6716/6731 line ratio is close 
to the limit for which this ratio is not a reliable indicator of the 
density. A more accurate estimation of in this knot should be 
derived from the [Fen] 7155/8617 line ratio. However, although 
the [Fen] /i 7155 A line is detected in knot A, unfortunately the 
[Fen] A 8617 A line is not included within the spectral range 
used. It should be remarked that such a high is most probably 
real (i. e. is not an artifact: first, because a similar low ratio has 
been found in each of the three adjacent pixels binned to get the 
A2 spectrum; and second, because the same behaviour at this 
knot position was obtained from another exploratory spectrum 
acquired one year before). Thus, this result gives also support 
to the idea that the knot A emission is split into two different 
substructures that are not well resolved by our data. 

A similar analysis, but for the slit position S2, intersecting 
emission from knots C to F, is shown in Fig.fTO] Regarding the 
gas excitation, the [N n]/Ha and [S ii]/ Ha line ratios suggest that 
three different regions should be defined along the slit: the first 
one, extending up to ~ + r.'5 from the knot C peak, having the 
highest [Sii]/Ha ratios (> 1) (the lowest gas excitation found 
along S2, being however higher than the excitation infered in 
knot B); the second one extending up to 7" west beyond the knot 
C intensity peak and enclosing knots D and E, where the exci- 
tation increases, although it is lower than in knot A; and finally, 
the western region beyond the end of knot E, which includes 
the F filament and the Fl faint knot, where the gas excitation 
decreases again. The spatial behaviour of the electron density 
(derived from the [Sn]A 6716/6731 ratio) along knots C to F 
is displayed in the bottom panel of the figure. The trend found 
is sugestive of density inhomogeneities in the emitting gas, the 
sizes of the clumps remaining unresolved with our spatial reso- 
lution. As in the case of the knots A and B emission, we found 
two densitiy gas components: a lower-density emitting gas, with 
100 < «e< 300 cm"^ and a higher-density component, with 700 
< 2500 cm"-^. A more careful inspection of the We spatial 
distribution also shows that there are several density enhance- 
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Fig. 8. Spectra at selected positions along the slits, from east (top left) to west (bottom right), belonging to the knots labeleld in each 
panel. Each spectrum was obtained by binning the signal within a window of three pixels (~ 0'.'6) centered at the position labeled in 
the panel. The reference postition (Off=0") for the slit S 1 corresponds to the [S ii] peak intensity position of knot A, which coincides 
with the substructure Al. The reference postition (Off=0") for the slit S2 corresponds to the [S ii] peak intensity position of knot C. 
In all the spectra, the intensity has been normalized to the Ha peak intensity. 



ments within each knot, possibly indicating that the knots could 
be split in substructures by imaging the emission with higher an- 
gular resolution. The [N ii]/[0 1] line ratio only depends weakly 
on the temperature, being a good indicator of the leve l of gl obal 
ionization in the flow (see e. g. Bacciotti & Eisloff'el [19991). Its 
values as a function of the slit position are plotted in Figs.|9]and 
[TOlto qualitatively show the variation of the degree of ionization 
in the outflow. As can be seen from these figures, this ratio is ap- 
preciably higher (> 1) along knot A than along knot B, (where 
the ratio is < 1, the lowest value found for all the knots). Thus, 
knot A appears as the most excited and highly ionized region 
in HH 223, while the closest region, knot B, is the lowest ion- 
ized a nd excitated one. From the Hartigan, Morse & Raymond 
( 1199 4^ models, the [N n]/[0 1] line ratios allowed us to constraint 
the ionization fraction (jCg) to 25-30% along knot A, 15% along 
knot B and 20-25% along knots C to E, all these Xg values be- 
ing roughly independent of the range of preshock magnetic field 
strengths expected within molecular clouds. The [N ii]/[0 1] ra- 
tios also constraint the shock velocities to 60-70 km s"'. From 
the estimated values and the derived n^, we also estimated the 
density (nH= njx^) for the knots. We found that knots A and 
B are significantly denser ( nH> 10"^ cm"^) than knots C to E, 
where hh ranges from IxlO-' to 3x10^ cm"^. 



3.2.2. Small-scale radial velocity along HH 223 

The radial velocities along HH 223, with a spatial sampling of 
~ 0'.'6, were derived from the line centroids of Gaussian fits to 
the Ha and the [S ii] A 6716,6731 A lines. Results are displayed 
in Figs. [TT| and [T2]for the slit positions SI and S2 respectively. 
From these figures, the following trends are found: i) The ve- 
locity behaviour as a fuction of the slit position is very similar 
for the lines considered, ii) The radial velocities derived appear 
highly blueshifted at all the positions enclosed within the knots 
A to E and Fl; velocities values slightly blueshifted or at the 
ambient gas velocity are found at positions corresponding to the 
F filament, but out of the Fl faint knot, iii) In several regions, 
significant variations in the velocity can be found at scales < 1", 
smaller than the knot sizes. 

Figure[TT|displays the radial velocity behaviour along knots 
A and B. There are appreciable velocity changes with position 
along knot A, the velocities being more blueshifted towards the 
east. Such a behaviour is found for all the emission lines (see 
also Fig. |5j. We derived a difference in velocity of ~ 60 km s"' 
between two position within knot A, which corresponds to the 
two emission peaks, Al and A2, of the narrow -band image (see 
also Fig. Q. The radial velocity has only small variations along 
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Fig. 9. Physical conditions along knots A and B (slit position 
SI). Upper panel: Contour plot of the Ho- emission, as in Fig.|7] 
The slit position has been marked by the horizontal lines. Lower 
panel, from top to bottom: excitation ([Nii]/Ha, [S ii]/Ha), ion- 
ization ([Nii]/[Oi]), and (in logarithmic scale, obtained from 
the [S ii] A 6716/6731 ratio), as a function of the position along 
the sUt. The reference position for the offsets is as in Fig.|5] 



knot B, reaching similar values to those found for positions east- 
wards of Al. 

Figure [12] displays the radial velocity behaviour along S2. 
The more blueshifted velocities (~ -130 km s"') are found for 
the knot D emission peak position, reaching similar values to 
those found for knot B. Less blueshifted velocities are found 
eastwards of knot D, reaching a value of ~ -60 km s 'around 
the knot C emission peak position. There is an appreciable 
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Fig. 10. Same as Fig.|9]but for the slit position S2, intersecting 
emission from knots C to F. The reference position for the offsets 
is as in Fig.|6] 



change in velocity for positions offset ~ +8" to +11" from knot 
C, where the derived velocities are compatible with the Vlsr 
value. Interestingly, the emission intersected by the slit at these 
positions mostly comes from the F filament, out of the three faint 
knots (F1-F3) found engulfed in it. Thus we speculate that the 
emission in this region probably would arise from ambient gas 
of a wall cavity which is already been excited and dragged by 
the supersonic outflow. 

Finally, it should be pointed out that, at several positions, 
the line profiles in all the strong emission lines of the unbinned 
spectra are suggestive of being double-peaked. Due to the spec- 
tral resolution used, we have not been able to resolve with con- 
fidence the two velocity contributions suggested by these line 
profiles. However, we believe that the contribution from two ve- 
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Fig. 11. Lowel panel: Radial velocity along knots A and B, de- 
rived from the Ha and the [S ii] A 6716, 6731 A line centroids. 
Upper panel: contours plot of the Ha image and the S 1 slit pro- 
jection. 



locity components is reliable and should be characterize from 
higher spectral resloution data. 



4. Conclusions 

In order to get knowledge on the kinematics and physical con- 
ditions in HH 223, we performed long-slit spectroscopy at two 
slit postions, covering emission from all the bright knots (A to 
E) and from the low-brightness F filament. 

The spectrum of each of the HH 223 bright knots was ob- 
tained by integrating the emission within the slit window aper- 
ture encompasing the spatial extent of the knot (ranging from ~ 
2" to 3'.'5). The analysis of these spectra allowed us to character- 
ize the nature of the knot emission. The results are summarized 
as follows: 

- For all the knots, the spectrum appears as produced by shock- 
excited gas, characteristic of the HH objects. 

- All the knot spectra have an intermediate/high degree of ex- 
citation, as their [Sii] and [Nii] to Ho- ratios indicate. The 
estimated ionization fraction (xg) ranges from 0. 15 in knot B 
to 0.3 in knot A. The electron density, derived from the [S ii] 
6716/6731 ratios, ranges from ^ 240 (in knot C) to 2800 
cm""*. The total density ranges from nu> lO"* cm"^ in knots 
A and B, to «h - lO'' cm"^ in knot E. Knot A is the dens- 
est and has the highest excitation and ionization. The lowest 
excitation and ionization is found in knot B. 

- The kinematics derived from the emission lines indicates that 
the knot emission is supersonic, with blueshifted velocities 
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Fig. 12. Same as Fig.[TT]but for the slit position S2, intersecting 
emission from knots C to F filament. 



ranging from -60 to -130 km s the more blueshifted ve- 
locity values being found in knot B. Velocity values compat- 
ible with the ambient gas are only derived for the F filament. 

From the spectra obtained by binning the signal of each three 
adjacent pixels (~ 0'.'6) along the slit, we searched for variations, 
both in the kinematics and in the physical conditions, at a scale 
smaller than the knot sizes. The relevant results are the follow- 
ing: 

- Reliable variations at a scale smaller than the knot sizes are 
found in both, kinematics and physical conditions. 

- Emission from knot A can be resolved into two substruc- 
tures (Al and A2), with a projected offset of ~ CK'S be- 
tween their Ha emission peaks, and a difference in velocity 
of ~ 60 km s"', the eastern (Al) substructure being the most 
blueshifted one. 

- Spatial variations of «e up to one order of magnitude are 
found. The lower «e values are found at positions where the 
slit mostly intersects gas from the low-brightness filamen- 
tary nebula (F). The highest «e is found coinciding with the 
western substructure (A2) of knot A. Other enhancements of 
«£ are found, in addition, at several positions inside knots C, 
D and E. 

In summary, concerning the nature of the emission, the 
spectroscopy let us to conclude that the emission from the 
HH 223 knots arises from shock-excited, highly blueshifted gas, 
being thus characteristic of an optical outflow having an in- 
homogeneous, knotty structure in which each knot could en- 
close smaller subcondensations, as suggested by the small-scale 
changes found in the kinematics and physical conditions. The 
knots emission should then trace internal working surfaces of 
shocks that originate because the supersonic gas is being ejected 
at varying speeds or with different ejection directions, and faster 
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ejecta overtake slower ones. A variable ejection velocity and a 
wiggling knot pattern is consistent with the exciting source being 
a YSO binary system. The H2 emission detected at the HH 223 
location should arise from gas collisionally excited by shocks, 
which is not dissociated because its lower velocity in the shock 
frame. Furthermore most of the emission from the HH 223 low- 
brightness optical nebula should arise from the gas excited and 
dragged by the CO outflow and traces an accelerated layer in the 
cavity walls, since velocities compatible with the ambient gas 
velocity are only derived at positions on the lower-brightness 
optical nebula (F filament), out of the Fl faint knot. 

In order to complete the knowledge of HH 223, it would be 
useful to map the whole extent of the emission with higher spa- 
tial and spectral resolution, searching for the knot substructures 
and resolving the contributions to the emission coming from the 
supersonic outflow and from the dragged, excited nebular gas. In 
addition, spectra through knot A, covering a wavelength range 
including the [Fe 11] 13F and 14F multiplets, are needed in order 
to properly derive rie around the A2 substructure from the [Fe 11] 
7155/8617 ratio, suitable to perform a «£ diagnostic when is 
higher than ~ 10'* cm"^, where the [Sii] 6716/6731 ratio is not 
density sensitive. 
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